Hello world

plot(runif(100))

dat<-read.csv("_dat/richness-data.csv")
head(dat)
##   X ele_meters richness   ln_ele ln_richness
## 1 1       1383 20.51757 7.232010    3.021282
## 2 2        764 15.61994 6.638568    2.748548
## 3 3        467 14.53794 6.146329    2.676762
## 4 4        457 10.15791 6.124683    2.318253
## 5 5       1779 17.97125 7.483807    2.888773
## 6 6       1415 13.06136 7.254885    2.569658
plot(richness~ele_meters,dat) 

fit<-lm(ln_richness~ln_ele,dat)
fit 
## 
## Call:
## lm(formula = ln_richness ~ ln_ele, data = dat)
## 
## Coefficients:
## (Intercept)       ln_ele  
##    -0.02814      0.40713